Spectrum and Dynamics of the BCS-BEC Crossover from a Few-Body Perspective. 
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The spectrum of two spin-up and two spin-down fermions in a trap is calculated using a correlated 
gaussian basis throughout the range of the BCS-BEC crossover. These accurate calculations provide 
a few-body solution to the crossover problem. This solution is used to study the time-evolution of 
the system as the scattering length is changed, mimicking experiments with Fermi gases near Fano- 
Feshbach resonances. The structure of avoiding crossings in the spectrum allow us to understand the 
dynamics of the system as a sequence of Landau-Zener transitions. Finally, we propose a ramping 
scheme to study atom-molecule coherence. 
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routinely explored in experiments with ultracold Fermi 
gases [3, 0, la 0, & El ■ I n this Letter, we explore the spec- 



Optical lattices are a powerful tool to study few body 
systems. When tunneling is negligible, optical lattices 
can be viewed as an ensemble of individual harmonic 
traps where the properties of these systems can be stud- 
ied. The interaction between the particles can be tuned 
using a Fano-Feshbach resonance [l| and the number of 
particles in each lattice site can be controlled 0, El • In 
a recent experiment with optical lattices, the spectrum 
of two fermions in a trap has been measured 0], demon- 
strating that few body trapped systems can be studied in 
their own right. Also, the BCS-BEC crossover has been 
nely explored in 

aaffaii. 

trum and dynamics of four trapped particles and we show 
how a few-body formulation allows us to obtain accu- 
rate solutions of the system without making the standard 
approximations of many-body theory. This provides an 
explicit representation of avoided crossings between the 
atomic degenerate Fermi gas (DFG or BCS)-type states 
and molecular BEC-type states. Our results directly ap- 
ply to optical lattice experiments and they provide a few- 
body perspective on BCS-BEC crossover dynamics. 

Specifically, we calculate the spectrum of two pairs of 
trapped fermionic atoms interacting through short-range 
potentials, all with the same mass m. One pair is as- 
sumed to be distinguishable from the other pair, but the 
two atoms within each pair are indistinguishable. The 
s-wave scattering length a of the short-range interactions 
will be tuned in the standard manner, [lj which allows us 
to explore the BCS-BEC crossover as a function of inter- 
action strength near a broad Fano-Feshbach resonance. 
Even though the BCS theory is not expected to apply 
to a 4-particle system, we still use this term to refer to 
the dynamical regime where a is small and negative. By 
solving the problem from a few-body perspective, we are 
able to give accurate properties - especially energy levels 
as well as time-dependent dynamics - of the full quantum 
mechanical spectrum at zero temperature. As a result we 
achieve a deeper understanding of the global topology of 
the spectrum, in addition to making quantitative predic- 
tions of transition probabilities and dynamical properties 
of this system when interactions change with time as in 
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To obtain the energ y sp ectrum, we use a correlated 
gaussian basis set [lfj, 111]. A diabatization procedure 
reduces the system to a tractable number of relevant 
eigenfunctions, after which we solve the time-dependent 
Schrodinger equation using the diabatic representation. 

The Hamiltonian adopted is 
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In this convention, particles 1 and 2 are "spin up" and 
particles 3 and 4 are "spin down" . The two-body po- 
tential function V(ry) is taken to be a purely attractive 
gaussian, V{vij) = Vq exp(— r?-/2dQ) where the width do 
of the gaussian is fixed and the depth Vq is tuned to 
produce the desired two-body scattering length a. To 
obtain results independent of the model potential prop- 
erties, we concentrate on the range do <C a; lo , where 
a ho = (h/mujo) 1 ^ 2 is the trap length. By considering 
different widths do ranging from 0.05a/! O to 0.01a^ o we 
have verified that our results exhibit a weak dependence 
on do- All results presented in this Letter correspond to 
d = 0.01a^ o . The eigenspectrum of Eq.fll]) is obtained 
by an expansion into correlated Gaussian basis set, i.e., 



^{ dij }(ri,r 2 ,r 3 ,r4) = S {Vo(Rc7M)e~ 



(2) 

with Rcm = (ri + *2 + i"3 + i"4)/4 the center-of- 
mass coordinate, ipo the center-of-mass ground state 
■)/>o(Rcm) = e~ 2RcM / ah ° , and S the symmetrization op- 
erator. Opposite-spin fermions are treated as distinguish- 
able particles, so S = (1 — 7 7 i,2)(l — ^3,4) where V is 
the permutation operator. The wavefunctions obtained 
from this basis set are in the ground center-of-mass state 
(with J cm — 0). The relative coordinate wavefunction 
has quantum numbers J T = + , since the basis functions 
only depend on Rcm and the interparticle distances, ry. 

The basis functions, defined in Eq.((2]), are character- 
ized by the set of values {<%} which are selected semi- 
randomly. While the dij corresponding to different spin 
fermions are selected to range from a fraction of do up to 
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a couple of times a^ to describe dimer formation, those 
dij corresponding to equal spin fermions are selected to 
be of the order of the trap length aho- The typical size 
of the basis set used in these calculations is about 7000- 
15000. The advantage of the correlated gaussian basis 
set is that all the matrix elements can be calculated ana- 
lytically in terms of the {dy}'s and the properties of the 
two-body potential and the trap. Prior to diagonalizing 
the Hamiltonian, a linear transformation eliminates lin- 
ear dependence, typically reducing the basis set size by 
less than 10%. The basis set is fixed while Vq is tuned 
to give different scattering lengths, whereby matrix ele- 
ments are calculated once and then used to obtain the 
spectrum throughout the entire range of the BCS-BEC 
crossover. 

The accuracy and convergence of the calculations have 
been verified in detail. The ground state energy agrees 
with fixed-node diffusion Monte Carlo (FN-DMC) cal- 
culations throughout the BCS-BEC crossover [12]. For 
example, the ground state energy at unitarity varies from 
5.099huj to 5.027huj as we change do from 0.05a/ lo to 
0.01a/j O ; FN-DMC calculations for a square well potential 
of range O.Ola/jo lead to a ground state energy 5.069(9)fiw 
fljj ]. The difference between our results and the FN- 
DMC results is thus about 1%, where the shape and range 
of the interactions start to play a role. Higher excited 
states are in agreement with the BCS and BEC limiting 
behaviors. Furthermore, the ground and excited states in 
the BEC can be used to extract the dimer-dimer scatter- 
ing length, which agrees with the Petrov et al. prediction 
[ 1 31 ] . and the corresponding effective range 12 1. 

The spectrum as a function of A = 1/a shows a series of 
apparent crossings and avoided crossings in the unitarity 
region. The avoided crossings can be roughly character- 
ized by their width AA, the range over which the two 
adiabatic eigenstates interact appreciably, into two main 
categories: narrow crossings, where AA <C l/a/ lo , and 
wide crossings, where AA > l/a^o- We adopt a variant 



of the diabatization procedure presented in Ref. [lj] to 
diabatize narrow crossings while leaving wide crossings 
adiabatic, which gives smooth energy curves. Figure [T] 
presents the partially-diabatic spectrum in the BCS-BEC 
crossover. The inset shows a zoom of the transition re- 
gion, i.e., the strongly interacting regime where the avoid- 
ing crossings occur. This structure of avoiding cross- 
ings permits a global view of the manner in which states 
evolve from weakly interacting fermions at a < to all 
the different configurations of a Fermi gas at a > 0, i.e. 
molecular bosonic states, fermionic states and molecular 
boson- Fermi mixture (see Fig. [I]). Furthermore, it allow 
us to visualize concretely the possible pathways of the 
time-dependent sweep experiments, as is shown below. 

A partially-diabatic representation can be used to de- 
scribe a ramp of an initial configuration through the 
BCS-BEC crossover as in the experiments carried out 
at different laboratories, like JILA and Rice. The ini- 




FIG. 1: (Color Online) The energy spectrum for four particles 
in a trap is shown versus the dimensionless quantity Xaho (20 
states that correlate diabatically with the 20 lowest energies in 
the noninteracting limit). The black solid curve is the ground 
state. The blue dashed curves are states that diabatically 
approach excited dimer-dimer configurations, the red dash- 
dotted curves correspond to states that represent one dimer 
plus two free atoms on the BEC side, while the states in 
green circles correlate diabatically to four free atoms. The 
lowest curve drawn with green circles is the Fermi gas ground 
state on the BEC side of the resonance. Inset: zoom of the 
adiabatic spectrum in the crossover regions. In this figure, all 
states are considered, showing the rich structure of avoiding 
crossings. 



tial configuration is propagated using the time dependent 
Schrodinger equation. Starting from the ground state in 
the BCS side, the parameter A is ramped through the res- 
onance to the BEC side at different speeds v = 4r . In a 
homogenous system, the ramps are only characterized by 
the initial density p, the speed v and to. This suggests 
that x = f^lisf I is the relevant dimensionless quantity 
that characterize the ramp speed in large systems. Thus, 
we will use x to compare 4-body results with large sys- 
tems. This idea of using the density, and not properties 
of the trap, to connect few-body calculations with large 
systems has been previously implemented [15 1. 

To interpret our numerical results we apply the 
Landau-Zener approximation, which predicts that the 
probability for a transition from the adiabatic tyj to ^ 
is Tij(x) — e~* ij / x , where Xij are dimensionless param- 



eters extracted from properties of the adiabatic eigen- 
states, as is discussed below. 

The nonadiabatic coupling or P-matrix controls the 
probability of nonadiabatic transitions to a good approx- 
imation using the Landau-Zener model. The coupling 
between two adiabatic states and $>j is Pij(X) = 
(^il^gjf-) where A is the adiabatic parameter. Clark 
has shown that, if the transition has a form consis- 
tent with the Landau-Zener approximation, then the P- 
matrix element for a transition from $j to 'J'i has a 
Lorentzian form whose width, along with the correspond- 
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ing eigenenergies, characterize the Landau-Zener param- 
eter Kij [l6| . We numerically evaluate all the potentially 
important and verify that the largest couplings gen- 
erally have a smooth single-peak form that is approxi- 
mately Lorentzian. The largest Pij relevant for this spe- 
cific dynamical sweep correspond to transitions among 
Wi, ^2, ^5 and ^14. Here, the partially-diabatic states 
are labeled in increasing order of their energy on the BCS 
side (see Fig. [I]), refers to the ground state, \I/ 2 is the 
first excited dimer-dimer configuration, \tg is the second 
configuration with a dimer and two atoms, and ^14 is 
the lowest configuration with no dimers, the "fermionic 
ground state on the BEC side" . The final probability dis- 
tribution can be explained as sequence of Landau-Zener 
transitions between these four partially-diabatic states 
where the positions of the peaks of the Pij (A) determine 
a specific order in which the transitions occur, namely 
1 — > 2, followed by 2 — ► 5 and 1 — ► 5, and finally 5 — > 14. 

We use pi to denote the probability of ending up in 
state following one ramp that started out in the DFG 
ground state (#1) on the BCS side. This Landau-Zener 
model then predicts that 

p l = (l-T B) i)(l-T 2 ,i), 

P2 = (1 - 75,2)T 2) 1, 

P5 = (1 - T 14>5 ) [T M (1 - T 2 ,i) + T 5j2 T 2>1 ] , (3) 

Pl4 = Tl4,5 [75,l(l -72,l) +75 ;2 T 2 ,i] . 

The sum of all these probabilities is unity by construc- 
tion. The Landau-Zener parameters obtained from the 
P-matrix analysis are >£ 2 ,i ~ 35, xs t i as 43, X5 i2 as 14 
and >£'i4 i 5 as 90. Interestingly, the Landau-Zener model 
shows decent agreement with the numerical results de- 
spite the neglect of many possible transitions. Figure [2] 
presents the probability of evolving in a certain configura- 
tion following a single unidirectional ramp, as a function 
of the dimensionless ramp speed parameter \. The nu- 
merical ramps are initiated at A^ ~ —7/aho and finalized 
at A/ ~ 7/aho- As the speed is increased, the final state 
changes from a molecular bosonic ground state, i.e., ^1, 
to a "fermionic ground state", i.e., ^14, in qualitative 
agreement with experiments. 

To relate our results with experiments carried out at 
JILA and Rice, we write the Landau-Zener parameter for 
atom-molecule transitions 5 — x mo i fx m terms of exper- 
imentally accessible variables. If the dependence of a on 
the magnetic field is approximated in the usual manner as 

a(B) = ab g ^1 + B ™ Bo ^ then v = (dB fdt)/wab g (mea- 
sured near unitarity, i.e., where the transitions occur). 
Therefore, S — x mo i{dB / dt)~ l ph\wab a \ /m which agrees 
with previous theoretical predictions 17, [ijj]. The de- 
pendence of S on p has been experimentally verified [l^ | . 
To evaluate x we use the average density for the non- 
interacting 4-particles Fermi gas, namely p as 0.153/a^ o . 
The molecular fraction, the fraction of atoms that be- 
come molecules after the ramp is over, is probably the 
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FIG. 2: (Color Online) The probability of evolving into a 
given configuration is shown as a function of the dimension- 
less ramp speed parameter \- The symbols correspond to 
the full numerical solution while the curves are Landau-Zener 
results. The black solid curve and circles correspond to the 
ground state configuration. The blue dash-dotted curve and 
crosses correspond to higher dimer-dimer excitations. The red 
dashed curve and squares correspond to ramps that produce 
a dimer plus two free atoms. The green dotted curve and 
diamonds correspond to the lowest configuration of four free 
atoms, i.e., the Fermi gas "ground state" on the BEC side of 
the resonance. 



most relevant quantity to compare with experiments. For 
our 4-body system, the molecular fraction is defined as 
the probability of ending up in a dimer-dimer configura- 
tion plus half of the probability to form a configuration of 
the "dimer plus two free atom" type. In the experiments 
we compare with, the molecule fraction was fitted to a 



Landau-Zener function, i.e. 
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where f m is is the maximum conversion efficiency which 
depends on temperature. 

Whether a Landau-Zener function is the correct func- 
tional form to describe the molecule formation fraction 
in large systems remains a question that existing experi- 



ments have not resolved 3 2(| 21 1. The Landau-Zener 
model presented in this work for four particles does not 
predict a single Landau-Zener function but a combina- 
tion of different Landau-Zener terms. However, the final 
molecule fraction predicted by this model and the nu- 
merical results for the molecule formation fraction can 
be approximately fitted by this Landau-Zener function 
59 ± 6; this value is higher than the two- 



with x' 



(4) 



body prediction of x. 



(2) 



42. Our results are consistent 



with the experimental Landau-Zener parameter obtained 
in Ref.0] for 40 K. The fit of the experimental data to a 
Landau-Zener formula predicted a as 62 ± 15. Also, 
experiments carried out at Rice measured the Landau- 
Zener parameter for 6 Li 7]. Taking into account the 
conditions of the experiment and the properties of the 
6 Li Fano-Feshbach resonance at B as 543. 8G, we esti- 
90. Both experiments were carried out at 



mate x. 



mol 



finite temperature and, consequently, maximum conver- 
sion efficiencies were approximately f m ~ 0.5, while our 
calculations are at T = where f m = 1. In summary, 
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both experiments are in general agreement with the four- 
body predictions. In the following, we propose more so- 
phisticated sweeps to study atom-molecule coherence. 
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At (l/w ) uj (lo ) 



that a simple Landau-Zener model approximately de- 
scribes the dynamics of unidirectional ramps. The spec- 
trum and dynamics of this system is interesting for op- 
tical lattice experiments. These would allow access to 
physics that cannot be probed in the two-body system, 
like corrections to the energy spectrum due to the dimer- 
dimer interaction, and also to the atom-dimer interac- 
tion. Also, the system of two "spin up" and two "spin 
down" fermions in a trap exhibits many of the ingredients 
of the BCS-BEC crossover problem, and in that sense the 
present results provide a few-body perspective on Fermi 
gas experiments. 



FIG. 3: (a) Molecular formation fraction is shown as a func- 
tion of the delay At. (b) Fourier transform of the Fig. 3(a). 
The peaks of the spectrum correspond to the most important 
energy transitions. 
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Atom- molecule oscillations [22( or quantum beats [15 1 
have previously been explored for condensates, and 
also for fermionic systems near a narrow Fano-Feshbach 
resonance (23j|. To study atom- molecule coherence in a 
Fermi gas near a broad Fano-Feshbach resonance, we 
have considered different ramping schemes and noticed 
one that enhances the atom-molecule oscillations. Start- 
ing in the ground state on the BCS side, one ramps at 
medium speed (x ~ 80) to the BEC side and pauses for a 
time At at a value X s uu ; then one ramps back at the same 
speed to the BCS side where the scattering length is close 
to zero. Finally, one slowly ramps A to the BEC (x ~ 7) 
side and measures the resulting molecular fraction. This 
is shown as a function of At in Fig. Ela). Observe that 
this ramping scheme produces large coherent oscillations 
in the molecular fraction. To interpret their frequencies 
we Fourier transform the time-dependent molecular frac- 
tion (see Fig. OUJb)). The frequency domain peaks corre- 
late with the most important configurations during the 
waiting period at \ s uu ~ 5/a/ lo . The Bohr frequencies at 
uj an 28wo correspond to coherences between states differ- 
ing in one broken dimer bond. For example, the highest 
peak is a coherence between VE^ and $14 while the second 
highest is a coherence between ^1 and $5. The frequen- 
cies around lo w 57luq correspond to coherences between 
states differing in two broken dimers bonds, e.g. coher- 
ences between vf^ and ^14, and between ^2 and ^14. In 
optical lattice experiments, this kind of multipeak struc- 
ture should be particularly pronounced in a tight trap. 

The four body problem remains fundamental and chal- 
lenging. We have presented an accurate numerical solu- 
tion of the spectrum of two "spin up" and two "spin 
down" fermions in a trap throughout the range of the 
BCS-BEC crossover. Even though the spectrum presents 
a rich structure of avoided crossings, we have shown 
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